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We have proved 


THEOREM 2. The rectangular distribution function of refractory 
periods 


N (6) = 


1 
fon 0i< 0.0905 
02 — 01 


= 0 elsewhere 


gives rise to an output intensity function which tends asymptotically 
to the value (62 — 6,)7 log (62/61) as x grows without bound. 

The asymptotic value of the output intensity determined by the 
distribution function (14) is, of course, a function of the parameters 
of the aggregate, namely, the minimum and the maximum refractory 
periods in it. If we hold the minimum fixed, we can see how the out- 
put intensity varies with the range (or spread) of the refractory 
periods in the aggregate. We see that as the range narrows down to 
zero, the limiting value of the output intensity becomes 


log —=— (19) 
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as should be expected, since giving every neuron in the aggregate the 
same refractory period is equivalent to considering the behavior of 
a single neuron. On the other hand, as the range increases, we have 


Lim 
6200 Oo a 3 


log — 0} dS 0: (20) 
L 01 

Thus the output intensity vanishes as the spread of refractory 
periods increases without bound. If, however, the thresholds of the 
neurons are proportional to their refractory periods, so that the 
total number of the neurons involved in responding to the input is 
proportional to the range of 6, then the total output F'(6,, 52), ac- 
cording to Definition 2 above, will be obtained by multiplying the 
output intensity by p(é6. — 6.) where p is a constant of proportion- 
ality. Then 


F (6;, 62) =p log =. (21) 
al 

This is a “mobilization” picture of growing excitation, where more 
neurons come into play as the input intensity increases. If 6, is fixed, 
and 6., being the upper limit of the refractory periods, proportional 
to the upper limit of the thresholds and therefore to the input inten- 
sity, then the total output as a function of input intensity follows the 
Weber-Fechner law. 

It is easy, of course, to get the logarithmic input-output curve 
in another, more obvious way, if one employs the pure “mobilization” 
model. 

In this model each neuron responds with a certain frequency 
characteristic for it independently of the input frequency, provided 
the latter exceeds a certain threshold. The increased total output of 
an aggregate accompanying an increased input is due entirely to 
“mobilization,” that is the involvement of more and more neurons 
as their threshold frequencies are exceeded by the input. Suppose 
now that the thresholds in the aggregate are so distributed that their 
frequency of occurrence is inversely proportional to their magnitude 
h in a finite range: 


N(h)=c/h for ky Shs he 
(22) 
= 0 elsewhere, 
where c = [log (h2/h1)]~ is the proper normalization constant. Then 
if the input intensity x is just sufficient to fire the neurons with 
threshold h(x), then all the neurons with thresholds below h will 
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respond with their characteristic output intensity K , which we here 
suppose the same for all neurons, and the total output of the aggre- 


gate will be 


» A(z) h 
F(a) =K | N(h)dh=A log —-, (23) 
7 1 
where A = Ke. If we further assume that h is proportional to x, 
(23) becomes 


F (x2) =A log (ka) (24) 


where i: is another constant. 

This model can be generalized by dropping the restrictions on 
h(x) and making the response frequency of each neuron a function 
of both its threshold and of the input frequency. In this case 


fal h(2) 
F(x) = | K(h,«)N(h)dh (25) 
hy 

which by a change of variables can be put into a form of Volterra’s 
integral equation of the first kind. Considerations involving refrac- 
tory periods can be incorporated into the kernel of the equation 
(which describes the behavior of the individual neuron). If the re- 
fractory periods are associated with the thresholds by some func- 
tional relation 6 = 6(h), we have an integral equation of an even 
more general type, 


S5[h(2)] h(@) 
F(x) = qf K'(h,6,2)N(h)N'(8)dhdd, — (26) 
5 (ha) Ia 
where N’(6) is the frequency distribution with respect to 6. Clear- 
ly, equation (2) is a very special case of (26). 
The Distribution N(6) = k? 6 e*, 


We shall examine one more type of distribution with respect to 
‘refractory periods, namely 


N(6) =k? 6 e™, (27) 
Here the output intensity is given by 


() 1 any 
f(e) =khe>d>— de*®dd, (28) 


n=1 
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which upon integration yields 
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If we let e*” =z, and thus x = — k/logz, the right side of (29) 
may be represented by 


(29) 


= z log(i — z) — log (1 — z) — log z log(1—z) 
z log z 


- 1 (30) 


as may be seen by expanding (30) in power series. We are interested 
in the limiting value of the expression (30) as « tends to infinity, 
that is, y tends to zero, and z tends to unity. After several appli- 
cations of |’HOpital’s rule, it turns out that the limit of the expres- 
sion in the bracket of (30) is —1, that is 
he 2 62 aw oe (31) 

Since the mode of 6 in the distribution (27) is re 1/k, we see that 
the entire aggregate behaves in the limit as a single neuron whose 
refractory period is 1/k. It can be shown, however, that the func- 
tion (30) is monotone increasing with respect to x , so that the oscil- 
lations which the output intensity of the single neuron exhibits are 
“smoothed out” in the aggregate. 

The same method can be used to compute the output intensity 
function and the limiting values thereof for other forms of the dis- 
tribution. 


Randomization of Stimulus Intervals. 


If we suppose that a neuron is a primary receptor, that is, re- 
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ceives stimuli directly from the outside, it is reasonable to random- 
ize the intervals between the incoming stimuli. For example the 
stimuli can be light quanta or some other form of radiation. 

Suppose that the distribution of the incoming stimuli in time is 
a Poisson distribution, that is, if an instant is picked at random the 
probability distribution of the time of the next stimulus will be given 
by 

P(2,7) =«xe", (32) 

where x is independent of time and ; is the interval between the 
chosen instant and the next stimulus. Clearly (32) is also the prob- 
ability distribution of the interval between two successive stimuli. 
Furthermore the expected interval is given by 


f are*tdr=1/z, (33) 
0 


so that x is seen to be the input intensity. 

R. Jost (1947) gives the relation between input and output in- 
tensities of a counter with a refractory period, which in our nota- 
tion is as follows: 


f(a) =a/ 1+ ip edt | (34) 


where “ is the average input intensity over the entire time range. 
If P(x, 7) is given by (32), © =a and (34) reduces to 


f(%) aos + 1). (35) 


However Jost’s formula holds in the general case where x is a func- 
tion of time. 

In the simple case of the Poisson distribution where x is inde- 
pendent of time, (35) can be derived from very elementary consid- 
erations. Note that the expected time of the next stimulus is 1/z 
regardless of what instant is chosen. In particular if the origin of 
time is chosen at the instant when the neuron fires and the expected 
time of the next impulse is measured from 6, the refractory period 
of the neuron, then the expected time of the next firing measured 
from the origin will be 6 + 1/z, and hence the output intensity will 
be given by (6 + 1/x)“, which is the right side of (35). 

It is important to note that the probability distribution of the 
intervals between successive firings will in this case not be of the 
form (82), since if a firing instant is chosen, the probability of the 
next firing will be zero for the next interval of time of duration 6. 
In fact in the extreme case, where the neuron is firing with its limit- 
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ing frequency 1/6, this distribution becomes the Dirac function with 
the peak at some point between zero and 6, depending on the instant 
chosen. If the instant chosen as origin is not specified (chosen at 
random), then the distribution is seen to be rectangular 


P(t) =1/6 for 0 St=6 

= 0 elsewhere. 36) 
The intermediate cases are represented by functions ranging from 
(82) to (86). 

We shall need these probability distributions in our discussion of 
filter nets. However their derivation will be postponed at this time, 
and we shall construct the simplest sort of a filter net, where the 
refractory period of each neuron is zero. 


The Filter Net. a 


A filter net is an arrangement of neurons in which certain neu- 
rons are able to respond to only a finite range of input intensity of 
outside stimuli. Such a net is shown in Figure 3. 


Ficure 3. A simple filter net. 


Neurons N, and N,z in Figure 8 receive randomized stimuli from 
the outside. Neuron N, excites N;, but N. inhibits Nz; provided N3; 
receives a stimulus from N, within an interval o before it receives 
a stimulus from N,, in which case the stimulus from N, fails to ex- 
cite N,;. In short N; fires at the instant ¢, if it receives a stimulus 
from N, at that instant (synaptic delay is not considered) and if it 
has not received a stimulus from N, at the time 7, such that 
t—coZrSt. . 

Now the probability of receiving a stimulus from N, in the in- 
finitesimal interval dt is xdt , where x is the input intensity from the 
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outside. Moreover, because of (32) the probability of not having 
received a stimulus from N., during the interval o is e*’. This is so 
because xe*?dt is the probability of the immediately preceding im- 
pulse from Nz occurring between o and o + dt before t. Now «dt is 
the probability of a stimulus at an arbitrary instant. Therefore e°? 
is the probability that no stimulus has occurred within oc. 
Hence the probability of firing of N; at a given instant is given 
by 
f(x ,o0)dt=xe*°dt (37) 


and hence xe’ represents the output intensity. It is seen to be a 
monotone decreasing function of oc, but it has a maximum with re- 
spect to x , namely at x = o, where the output intensity is (ec)7. 


f W=x e7*% 


X 


1 
X; nr 
Figure 4. Neuron N, in Figure 3 responds only to outside input frequencies 
between x, and a, if 0 < h S (es)-1. 


Let now neuron N, have a threshold h, such that 0 < h < (@c)71. 
As can be readily seen from Figure 4, N, will respond to outside input 
frequencies only between the limits shown. As h increases and ap- 
proaches (éc)-* the range of response of N, diminishes, and the neu- 
ron gradually turns into a “discrimination” neuron, able to respond 
to practically a single frequency. If h increases still further, N, 
fails to respond altogether and the filter net breaks down. 

This work was aided by grants from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of The University of Chicago, and 


from The Lucius N. Littauer Foundation to The University of Chi- 
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It is shown that Blair’s theory of excitation is independent of, and 
consequently valid for, any possible relationship between the threshold 
and the magnitude of the stimulus. It is pointed out that if a depend- 
ence of threshold on stimulus is assumed, the concept of rheobase be- 
comes meaningless. Consequently, contrary to Blair’s impression, the 
disagreement between his theory in its original form and the experi- 
mental data on the time of incipient excitation with constant stimulus 
(response time) cannot be explained by assuming a dependence of the 
threshold on the magnitude of the stimulus. It is shown that a modifica- 
tion of Blair’s interpretation, obtained by taking into account effects 
of internal energy sources released by the stimulus, eliminates the dis- 
agreement mentioned above between theory and experiment. The role 
of such modification in connection with propagation of excitation is dis- 
cussed. 


The Blair theory of excitation is based, as it is well known, on 
the differential equation 


de/dt = ol — pe, (1) 


where J is the magnitude of the external stimulus (e.g. intensity of 
a constant current applied to a nerve), c is the magnitude of some 
effect of the stimulus in the biological system considered and o and p 
‘are two constants, the first measuring the susceptibility of the system 
to the stimulus J, the second its ability to react to the effect ¢ of the 
stimulus. If the magnitude of the stimulus J is constant, the explicit 
solution ¢(t) of equation (1), under the condition ¢(0) = 0, is: 


c(t) = (o/p)I(i—e*). : (2) 


The condition c(0) = 0 is inherent to the definition of ¢ itself. The 
assumption that a threshold magnitude h of ¢ is necessary to cause 
excitation leads to the following theoretical expression of the mini- 
mum value R of I required for excitation: 


R= (p/o)h. ? (3) 
This is the threshold stimulus, the so-called rheobase, whose appli- 
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cation for an infinite length of time would excite the system by 
bringing c from zero up to its threshold value h. If a constant stimu- 
lus I is applied which is larger than Ff, the duration 7 of applica- 
tion of such stimulus required to produce excitation is, from previ- 
ous equations: 

pt=In{I/U—f)]. (4) 


In the derivation of all the above relations the only assumptions 
made concerning the threshold h of ¢ are that this threshold exists 
and that it does not change during the application of the stimulus. 
No assumption of any kind is made as to a possible dependence of the 
threshold h on the magnitude I of the stimulus. Consequently, all 
the relations written above are valid without any change, even if h 
is assumed to be dependent on I in any fashion whatsoever. In other 
words, relation (4), exactly as it is written, is an unavoidable con- 
sequence of the theory. It is known that this relation does not agree 
entirely with experiments, which can be fitted, however, by the for- 
mula: 

pt—InfT/(I—R)] +C, (5) 


where C is a constant. Blair thought that this fact did not contra- 
dict his theory, but that it simply meant that the threshold h de- 
pended on the magnitude J of the stimulus. The argument which we 
just presented shows that his reasoning must contain an error. Ap- 
parently he did not notice that his general theory, including formula 
(4), is independent of, and consequently valid for, any possible re- 
lationship between the magnitudes of the threshold h and of the 
stimulus I. None of his formulae can be changed by assuming a par- 
ticular dependence of h on I, consequently none of such possible de- 
pendences can make an additive constant appear on the right-hand 
side of formula (4). 

The question arises—which exactly is the erroneous step in 
Blair’s reasoning? The answer can be given in the following form: 
If the threshold h of ¢ is assumed to be dependent on the magnitude 
I of the stimulus, the concept of rheobase, and consequently also 
any formula containing it, becomes meaningless. In fact, the rheo- 
base FR is defined as the minimum magnitude of a time-independent 
stimulus I causing excitation, ie. as the minimum value of I, con- 
stant in time, which is able to bring ¢ from zero to its threshold val- 
ue h. Therefore, if h is dependent on I, rheobase R is also depend- 
ent on I, that is, R is a function of I and consequently cannot be a 
minimum of I! 


Another form of presenting the argument is the following: Since 
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relation (5) is experimentally correct, one may ask what are the 
implications of formula (5) and of all the other formulae of Blair’s 
theory with exception of formula (4) ; in other words, what. are the 
consequences of equations (1) to (3) combined with equation (5)? 
To answer this we note that from formula (5) the following rela- 
tion can be easily derived by elementary algebra: 


ert" {1 — (R/1)]= (6) 
But from equation (2) we have by definition of h and +: 
h= (o/p)I(1—e*"), 
which becomes by equation (6): 
h= (o/p)I(1—e°) + (0/p) Re“, 
and by equation (3) 
h= (o/p)I(1—e*) +he*. 


From here h = (c/p)I, which, compared with equation (3), gives 
R =I, that is, any magnitude of applied stimulus would be the 
rheobase, which is impossible, because J is a variable and R is a 
constant by their physical and physiological meanings. This conclu- 
sion is logically inconsistent because it is derived from mutually in- 
compatible facts—equation (5) and Blair’s theory. It can be added 
to this argument that since Blair’s reasoning implies J = Rk, equa- 
tion (5) contains a term In o and is consequently meaningless. How- 
ever it should not be inferred from here that the error under dis- 
cussion lies in an inappropriate application of a mathematical pro- 
cedure; it is of a purely logical nature. 

While pointing out that equation (5) cannot be put into agree- 
ment with Blair’s theory by assuming a dependence of the thresh- 
old of ¢ on the magnitude of the applied stimulus, we should keep in 
mind that Blair’s theory describes a surprising amount of facts with 
extremely simple means. In view of our very incomplete knowledge 
of the nature of neurophysiological phenomena it is not surprising 
that the few elementary thoughts from which Blair succeeded in 
developing a theory make the latter sometimes diverge from experi- 
mental facts. Although he himself indicated at one time that the 
constant C may not represent any essential part of the mechanism 
of the process, but may be perhaps only an expression of a perturba- 
tion due to electrodes, it seems worthwhile to examine whether the 
disagreement under discussion cannot be connected with some more 
fundamental fact since many experiments have been analyzed in the 


126 ON BLAIR’S THEORY OF EXCITATION 


literature in connection with that constant C. 

Blair himself pointed out a possibility for modifying his theory 
so as to lead to relation (5) with C #0. In fact, if we assume that 
there is a latency or a time lag between the effect and the stimulus, 
we have to use a condition of the type c(t) =0 for t < t instead 
of the condition c(0) = 0 where ¢t,, the time lag, is a non-negative 
quantity. We then obtain from the differential equation (1), under 
the assumption of a constant stimulus 7, 


o(t) = (e/p)I[l— er? © y. (7) 
This equation, with the additional condition, 
c(t) —=0 “tor f=%, 


would now replace relation (2). The rheobase would be still given 
by equation (3) which follows also from relation (7). But instead 
of equation (4) we would have: 


pr=It [1/7 — #)] + pt: 
A comparison of this with equation (5) gives 
C = P (i , 


which is positive. But in 54 out of 65 sets of experimental data ana- 
lyzed by Blair the constant C comes out negative. [The constant C 
in Blair’s papers equals —C of this paper. N. Rashevsky (1948, p. 
282) uses the symbol C, for our C; there is apparently a misprint on 
that page since the constant C, is negative for the data given on the 
same page.] Thus the concept of latency cannot be accepted within 
the framework of Blair’s theory. 

The constant C is obtained as an intersection of a straight line 
obtained by plotting experimental data on In [J/(J — R)] against 
7. Consequently experimental data on small values of 7 are impor- 
tant for the determination of C. As a matter of fact Blair’s data 
go much below 1 millisecond and the question arises whether dur- 
ing such short intervals of time after the closure of the circuit the 
transient may be sufficiently important so as to cause an appreci- 
able deviation of the intensity of the current from its asymptoti- 
cally constant value, since constant current is assumed in the the- 
ory. According to H. J. Curtis and K. S. Cole (1944) the membrane 
resistance of a squid giant axon is 1000 ohm/cm?, the membrane in- 
ductance is 0.2 henry/em®. The transient in an inductive circuit is 
given by: 
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[=[,(1-—e*'/*), (8) 


where J, is the asymptotic value of the current intensity on which 
theoretical calculations are based and 7 and L are the resistance and 
the inductance of the circuit. Assuming the above values of 7 and L 
we obtain from (8), for t = 4 millisecond: J ~ 0.918 J,. Thus the 
effect of the transient would amount to some 8% of the intensity of 
the current, if one neglects the effect of the part of the circuit ex- 
ternal to the biological system under experimentation. Since accord- 
ing to Cole et al (cf. Cole and Baker, 1941, p. 799; Curtis and Cole, 
1944, p. 805) the axoplasm and the intercellular fluid contribute to 
the resistance of the nerve but not to its inductance, it is likely that 
the actual transient is smaller than the above value. In this very 
crude estimate of the transient we have also neglected a possible 
capacity effect. It is easy, however, to show that a transient due 
to a small inductance would not lead to a formula of type (5) with a 
negative constant C. In fact, if we solve the differential equation 
(1) of Blair’s theory on the basis of the current (8) we obtain: 


1 (Innes E l I 
maul) te") 

Ras pike Pre disp 
where 1— r/L. Since 4 is very large, because the inductance is very 


small with respect to the resistance, we neglect the second term in 
the first parenthesis and obtain: 


ol, 
Cz = 
p 


The rheobase considered as the minimum value of J, causing excita- 
tion is still given by R = (p/c)h. Consequently, from equation (9) 
the following formula is obtained for the response time: 

i; A 


+ In —— 
I,—R wie. 


et | ; (9) 


A— p 


p7t=I1n 


with C = In [A/(A — p)], but this constant C is never real and nega- 
tive. 

The concept that energy stored in some form in the nerve is re- 
leased by the stimulus to work toward excitation is well known (Cur- 
tis and Cole, 1944, pp. 798-808; Lorente de No, 1947, pp. 440-441). 
It can be easily seen that such energy sources if taken into account in 
a suitable form together with Blair’s theory lead to equation (5) 
with a negative constant C as required by experiments. In fact, it 
is intuitive to assume that the amount of such internal energy re- 
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leased by the stimulus is proportional to the magnitude of the latter, 
within certain limits at least. Since Blair’s theory assumes that ap- 
plied current produces some effect in the nerve which leads to exci- 
tation, it is natural to assume that those internal energy sources 
produce the same effect when released, to some extent at least. We 
are lead in this way instead of equation (2) to the following rela- 
tion: 

c(t) = (o/p)I(1—e?' +a), (10) 
where a is a constant. The quantity (o/p)Ja measures the con- 
tribution of the internal energy sources to the effect c. The relation 
between the rheobase and the threshold is now: 


h= (o/p)R(1 + a), (11) 
which substituted into equation (10) gives for the response time 7: 
p7in f/U—R#)]—in(1 +a). (12) 


This is exactly equation (5) with a negative C, because a is posi- 
tive. In connection with equation (10) we have to assume that the 
internal energy is released during a certain very short interval of 
time between the moment of application of the current (¢ = 0) and 
the smallest response time 7 encountered in the experiments under 
consideration. Thus the release of energy whose effect appears in 
equation (10) would occur within a time tz smaller than 1/2 milli- 
second. Equation (10) would be valid for ¢ > t;. Since for the pur- 
pose of the present discussion we are not interested in magnitudes 
of c for t < tz, we can leave out of consideration the question of 
how fast that energy is released within that time t,;. 

To test this idea of internal energy sources in the form here pre- 
sented, let us examine whether it can be applied to the phenomena 
of propagation of excitation in a nerve fiber. We will follow the 
method applied by N. Rashevsky in chapter xxvii of his book (1948), 
in which he schematized the nerve fiber as a semi-infinite core con- 
ductor of radius S surrounded by a thin shell of thickness 6 and 
imbedded in a medium whose resistance is y times smaller than the 
resistance of the core of the same length. Under the assumption of 
a steady state, neglecting any possible contribution of internal ener- 
gies, the intensity of the current is given by the equation 


I= I, e**, (13) 


where I, is the intensity of the current at x =0 and 


honeets 1 
B are, ‘ 
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with r, the resistivity of the core and ry the resistivity of the shell. 
The condition of excitation is ] > R. Taking into account equations 
(11) and (18) it is seen that the excitation would occur only within 
the range of 2: 


ae oes Io o(1-+ a) 
p hp 

We can estimate this range by using the values of the constants based 
on H. A. Blair’s experiments as suggested by N. Rashevsky (1948, p. 
329),8 = 1,hp/ (Uc) =1/2. If there were no internal energy sources 
(a = 0) the excitation would propagate only within « < In 2 = 0.7 
em. A larger intensity of current Jp would be necessary to cause 
excitation beyond that range of x. With internal energy sources a 
is positive and the range is larger. For instance, with a value of 
C = —In(1 + a) = — 0.4, which is about the maximum value of 
C derived by Blair from his previously mentioned experiments, we 
obtain a = 1/2 and the range within which the excitation would pro- 
pagate comes out to « < 1.1 cm. It should be noted that a value of 
PB = 1 on which the above calculations are based implies a very low 
ratio of r./rs (cf. Rashevsky, 1948, p. 329). If this ratio is increased, 
6 is also increased and « is decreased. At the same time it should be 
kept in mind that a value of a =1/2, which represents the contribu- 
tion of the internal energy sources to the magnitude of c, amounts 
to 50% of the value that c would acquire asymptotically at t = 
if those energy sources were not acting [cf. equation (10)]. Thus 
there seems to be no reason to consider a higher value of a. Besides 
this, if the rheobase is determined directly by experiments the anal- 
ysis becomes completely independent of a. In fact, we obtain from 
equations (11) and (14) for the excited length of the nerve 


a= (1/6)In (1o/K). (15) 


It is seen in this way that taking into account the internal energies 
through a term a in equation (10), in agreement with some experi- 
mental facts, is however insufficient to explain the phenomena of 
propagation of excitation. 

The contribution of internal energies to propagation of excita- 
tion in the nerve has been taken into account successfully by N. Ra- 
shevsky in a different manner. We would like to show how his re- 
sults can be easily derived from the previous formulae. For this 
purpose let us note that by substituting equation (13) into (12) we 
obtain: 


(14) 


pr=—In(Ip—R ef?) + In, —In(1 + a) 
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which, differentiated with respect to 7, gives the velocity of propa- 
gation of excitation: 


& = (p/B) ((o/R) e®* — 1). (16) 


It may be easily seen from equation (15) that at the end of the 


range within which the excitation propagates the velocity x becomes 
zero as it should be, mathematically. We have seen that the shortness 
of this range contradicts the experimental facts. Rashevsky, however, 
takes into account the contribution of the internal energies by as- 
suming that the curve representing the intensity of the current, in- 
stead of being expressed at each moment by equation (13), is rigidly 
displaced along the nerve fiber with x = 0 always at the point of 
excitation. Under this assumption it is easy to see that excitation 
propagates with a constant velocity. In fact, if at a certain moment 
t the excitation reaches a point x, the velocity of propagation of ex- 
citation at that moment depends on the magnitude that ¢ has in an 
infinitesimal region dx immediately ahead of the excited point and 
on the magnitude that ¢ will have in the same region during an in- 
finitesimal interval of time dt immediately following ¢. But in that 
region and during that interval of time ¢ varies between [c(%, €) Jeze 
and 


, 0c(a, t) 2 c(x,¢) 
[ow t) owe + | ie dx + [9 i- dt, (17) 


where the subscript exc indicates that the functions in the brackets 
are calculated at that time ¢t at which the point 2 becomes excited. 
We will show that all the above terms in brackets are constant, which 
will prove that immediately ahead of the excited point, in time and 
in space, c has always the same value. This will prove that the ve- 
locity of excitation is a constant, because this velocity depends only 
on the magnitude of ¢c immediately ahead of the excited point. Now, 
by the definition of the process of excitation itself we have: 


[ec (x , t) ] exc — h= constant. 


To prove that the two derivatives in equation (17) are also constant, 
it is convenient to introduce a symbol for the distance covered by the 
excitation during a time ¢, viz.: &(t). The symbol & concerns in a 
sense the same quantity as x, except that x represents any point 
of the axis whereas &(t) expresses the relationship between that 
point and the time of its excitation. Thus x is an independent vari- 
able and €(¢) is a function. Since we disregard at present a possible 
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contribution of the internal energy sources through the factor a, we 
obtain from equations (12) and (13), putting x =a=0 : 


pt—=In [Io/ (fo — R)] ’ 


for the time of excitation t of the point x =0. Thus é(t) =0 for 
t < t. Assuming Rashevsky’s hypothesis of the rigid displacement 
of the current (13) together with the excitation, the explicit expres- 
sion for the intensity of the current as the function of x and ¢ be- 
comes: 

I= Ip e-B2+86(), 


Substituting this into the general formula for the solution of the dif- 
ferential equation (1), i.e. in: 


t 
c=ce?t i I(t’) et'dt’, 
0 


we obtain: 


0 


ot 
(at) =Ioe ere | ept’sBe(e”) Ge", (18) 


From here we easily derive 


dce(z,t)/ox=fe(z,t); 
d¢(x,t)/ot=—pe(x,t) + 1p o efé-02, 


But at the moment of excitation of x we have: 
c(z,t) =h, E(t) =z, 


the latter by the definition of the function é(t) itself. Thus the terms 
in brackets in equation (17) are constant in time and space and if 
we think of the excitation as a pulse propagating along the nerve, 
the situation ahead of the pulse would be always the same as it was 
at t = t, ahead of « = 0. Consequently the velocity of the pulse 
would be always equal to its initial value at x = 0. But this value 
is immediately obtained from the previous theory of excitation be- 
cause at the moment of excitation of x = 0 the quantity c is the same 
throughout the whole nerve fiber in both theories. We get in this 
way, putting « = 0 in equation (16): 


az = (p/B)[(Uo/R) —1]. 


This is Rashevsky’s formula. It may be added that this formula is 
in quite good agreement with experiments. In fact, taking into ac- 
count the expression of f, it is easy to see that the above formula 
implies a direct proportionality of the velocity of propagation to the 
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square root of the diameter. R. J. Pumphrey and J. Z. Young (1938) 
experimenting on nerve fibers of Cephalopods in the range of 30 to 
718u found a proportionality of the velocity to the power 0.614 of 
the diameter. J. B. Hursh (1939) and H. S. Gasser and H. Grund- 
fest (1939), experimenting with nerve fibers of cat and kitten, and 
cat and rabbit respectively, within ranges up to about 20u of diameter, 
found a proportionality of the velocity to the first power of the di- 
ameter or slightly less than that. The older experiments by H. S. 
Gasser and J. Erlanger (1927) imply also the same type of rela- 
tionship. The hyperbolic relationship between the velocity and the 
rheobase or the threshold is also confirmed by experiments (Fulton, 
1947, p. 102). Thus there are several experimental results which 
agree with Rashevsky’s formula. 

This work was aided by grants from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of The University of Chicago and 
from the William T. Morris Foundation to The University of Chicago. 
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Some consequences of Onsager’s generalization of Fick’s law are 
examined. It is found that metabolized solutes may flow continually 
against their concentration gradients. Inert solutes may exist in a 
higher or lower concentration inside of the cell than in the medium 
thus appearing to be accumulated or excluded. The magnitudes of such 
concentration differences are dependent upon the rates of metabolism of 
metabolized solutes. Alteration of these rates may further increase the 
concentration disparity by causing inert solute to flow from a low toa 
high concentration. Due to the mutual dependence, demanded by Onsag- 
er’s law, of the diffusion currents, the rates of chemically independent re- 
actions are mutually dependent. This so called coupling by diffusion 
implies that: The rate of metabolism of a given substrate is influenced 
by the rates of metabolism of metabolically unrelated substrates. Fur- 
thermore, the presence, in the medium, of an inert solute to which the 
membrane is permeable will influence the rates of concentration—de- 
pendent reactions in the cell. The spatial distribution of a catalyst in 
the diffusion field within the cell is examined. The general effect of 
including heat flow and thermal] diffusion in the cellular diffusion prob- 
lems is briefly pointed out. 


Introduction. It is well known that when two or more irrevers- 
ible transport processes (e.g. heat conduction, electrical conduction, 
diffusion) take place in a given system there is mutual interaction 
among the individual transports. Not only do the individual pro- 
cesses interfere with one another but in general they may not occur 
separately; i.e. one such process obligatorily entails the others. Vari- 
ous combinations of such “coupled” transports constitute familiar 
phenomena. 

An imposed electromotive force in a conductor of dissimilar 
metals causes the evolution (or absorption) of heat at the junctions 
(Peltier effect) ; the converse is the “thermoelectric force.” An im- 
posed electromotive force in an electrolytic solution causes the trans- 
port of matter; the converse is the electromotive force of a concen- 
tration cell. A solution in which there exists a gradient in tempera- 
ture exhibits also a gradient in composition. This phenomenon first 
studied by Soret is known as the “Soret effect.” The flow of matter 
due to a thermal gradient is called thermal diffusion; the converse 
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effect is less well known but has been qualitatively demonstrated 
(references cited in Onsager, 1931a). 

The above are examples of the mutual dependence of the trans- 
port of heat and electricity, matter and electricity, heat and matter. 
Usually the individual “flows” or “currents” are taken, in the ab- 
sence of interaction, to be given by empirical laws; viz., Fourier’s 
Law, Uhm’s Law and Fick’s Law. When interaction is not to be neg- 
lected the currents are taken to be given by phenomenological rela- 
tions which express any current as a linear combination of all of 


the gradients. For example for the current of diffusion, J, and of 


heat flow, q, the equations 
d=— (DV Gai Dav Ts 


GS (kVT ar keVC) 


express the assumptions that a temperature gradient produces a flow 
of matter and that a concentration gradient produces a flow of heat. 
The coefficients D , Dy; and k are the diffusion, thermal diffusion and 
thermal conductivity coefficients respectively; C is the concentration 
and T the absolute temperature. Relations analogous to the above 
but expressed in terms of the “forces” which produce the currents, 
viz. Vu and V(1/T) where uw is the chemical potential, are more 
often used (Onsager, 193la & b; Eckart, 1940; Leaf, 1946). When 
one considers multi-constituent solutions, there must be added to the 
above examples the mutual interaction among the transports of so- 
lutes (Onsager, 1931a). 


Onsager’s Generalization of Fick’s Law. 


The chief effect of the forces of interaction is of a thermody- 
namic nature: The force of diffusion, —Vyw, is different (usually 
less) from that which would obtain, under the same concentration 
gradient, in an ideal solution. Also the velocity of diffusion in a non- 
ideal solution is, due to electrostatic and “hydrodynamic” interac- 
tion, for a given force —Vu different from that in an ideal solution. 
The generalization of Fick’s Law (Onsager and Fuoss, 1932; On- 
sager, 1945) for a solution of s solutes relates the velocities of dif- 


fusion v; (relative to the solvent) linearly to the “forces” or poten- 
tial* gradients —V u;: 


*When electrolytes are invclved the u; are the electrochemical potentials, 
otherwise “; denotes the chemical potential, equation (6), of Gibbs. 
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Je UjC;—— > 0 Jigs @— 1, 2s (1) 


The solvent is assumed to be the dominant constituent; this is cer- 
tainly unambiguous when dilute aqueous solutions are being consid- 
ered, The flow of any solute is defined relative to a local coordinate 
frame moving with the solvent. Bulk motion (cf. also Eckart, 1940) 
determined by hydrodynamic considerations may also be imposed. 
One may in practice assume either that the velocity of the solvent 
or that the “bulk velocity” (Onsager and Fuoss, 1932), 


> ~ 


ISSI, Ve, (2) 


where V; is the partial molal volume of the ith solute, obeys the hy- 
drodynamic equation. The omission of the thermal diffusion term 
from equation (1) amounts to the assumption of an isothermal sys- 
tem in which there are gradients in composition and requires some 
comment, but this is deferred until later. 

The relations Q;; = 2;; hold for all 7 and 2; they state that the 
flow of solute 7 caused by a unit force on solute 7 is equal to the flow 
of 7 due to a unit force on 7. These symmetry relations, which seem 
physically necessary for the interaction expressed by (1), are im- 
plicit in the “principle of least dissipation of energy” (Onsager, 
1931a & b; Onsager and Fuoss, 1932) and are derivable on the basis 
of the principle of microscopic reversibility (Onsager, 193la & b; 
Casimir, 1945). 

The form of the simple Fick’s Law for a two component solu- 
tion analogous to equation (1) is 


es Cae (3) 


from which it is seen that 2 has the physical significance that Q/C 
is the “mobility” or the velocity caused by a unit force —Vuyu; the 
quantity C/Q is the “frictional coefficient” [cf. Hearon, 1950b, equa- 
tion (39), etc. and equation (8) below]. Comparison of (3) with the 
usual form 


jJ==—DVC (4) 
shows that 


D=2—. (5) 
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If wis given by 


= oral eee (6) 
where y is the activity coefficient, equation (5) gives 
olny ) 
= Ls Q (7) 
DG ind ( olnC 
For an ideal solution (y =1), (7) becomes 
O/C =O/h ie (8) 


It is of some interest to compare the remarks following (3) and 
equation (8) to the results obtained by G. Young (1938, pp. 171- 
172) from the Stokes-Navier equations and to those of N. Rashevsky 
(1949). 

In an ideal solution in which the solutes diffuse independently 
equation (1) reduces to 


Jp — Qa le (9) 


If departure from ideality is due to compound-formation, one dif- 
fusing substance carries another with it thereby influencing its flow 
and there will be coefficients Q;;, 7 # 7, which are not zero. On the 
basis of this physical picture the simple law (9) would hold for each 
compound and molecular species actually present in the system 
(Lamm, 1947). It is not always possible to satisfactorily describe 
departure from ideality in terms of compound-formation and the 
phenomenological description (1), albeit non-committal, is employed. 
Alternatively, if possible, the forces of interaction may be specified. 
Greatest progress along these lines has been made for the case of 
electrolytes, and Onsager and Fuoss (1932, p. 2761) give expres- 
sions for the order of magnitudes of the 2;;. Expressions for the 
forces of interaction and “drag forces” are also available for the im- 
perfect gas (cf. Young, 1938 and references therein; Landahl, 1942). 

While the form (1) is perhaps most useful for certain formal 
discussion we will find the alternative form, 


te Tie CN k=1,2,-:--s, (10) 


given by Onsager (1945), of greatest utility. He has presented re- 
lations among the D;; derived from equating the bulk velocity (2) 
to zero. A discussion of equation (10) similar to that presented for 
equation (1) can be given; relations between the Q,; and the D,; can 
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be derived (Eckart, 1940; Onsager, 1945). We will refer to the Des 
j #k, as the “drag coefficients” or “cross coefficients.” 


Purpose of the Paper. 
It is the purpose of this paper to present the solution of several 


problems in which the currents, J, , are taken from equation (10). 
Since some of these problems have been previously solved (Rashev- 
sky, 1948) under the assumption of the simple Fick’s Law, equation 
(4), it will be the purpose here to point out differences and new fea- 
tures introduced by retaining the cross-coefficients D;;, k # 7. It 
will be seen that the formal boundary value problem presented by 
even the simplest case is considerably different when the generalized 
Fick’s Law is used. Further, it becomes necessary to introduce gen- 
eralized permeability expressions and cross-permeability coefficients, 
hx; , analogous to the D,;. 

The situation arises as a consequence of (10), wherein it is pos- 
sible to have a stationary, non-equilibrium distribution of an inert 
solute in which the concentration is greater (or less) in the cell than 
in the environment. The disparity between the cellular and environ- 
mental concentration is determined by the rates of metabolism of 
whatever metabolites are being produced or consumed in the cell; a 
change in the rates of metabolism can cause inert solute to flow 
against a concentration gradient. This has been discussed briefly 
elsewhere (Hearon, 1950b) from the standpoint of thermodynamics, 
and the details of the diffusion problem are herein supplied. 

It can be seen readily from the results to follow that a conse- 
quence of equation (10) is, since the spatial concentration distribu- 
tion of a given solute is determined by the distribution of all other 
solutes, that the rate of one chemical reaction may be influenced by 
another of which it is quite independent chemically. This mutual de- 
pendence of chemically independent rates is called “coupling by dif- 
fusion.” 

Since in general the solutions of the diffusion equation are tract- 
able only when certain spatial symmetry can be assumed and, more 
especially, when stringently simplifying expressions are employed 
for the chemical rates of production and consumption, no very gen- 
eral treatment can be given. Considerable abstraction is necessary 
and, in fact, some physical approximation must be made. However, 
the salient features can be seen from the highly simplified (not very 
realistic) cases treated and certain qualitative generalizations are 
apparent. There are special and limiting cases of the problems dis- 
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cussed which are not investigated in detail but which are of interest 
for further study. 

The influence of the cross-coefficients upon the transient or time- 
dependent solutions is perhaps more interesting, and from the ex- 
perimental point of view more important, than the steady state situa- 
tions discussed here. Only the most simple variable state cases are 
briefly discussed in a subsequent note. 


The equation of continuity. 


By the superscript i we will denote quantities which refer to the 
interior of the cell and by the superscript e quantities referable to 
the environment. The equation of continuity for the kth substance 
inside of the cell is, in the steady state, 


—V-J> +q,=—0, kK=1,2,::::8, (11) 


where the J;' are obtained from (10) by inserting the values D,,;' 
and C;’. Here the first physical approximation appears: We are 
neglecting thermal diffusion and the “inverse Soret effect”? and as- 
suming an isothermal system. This is discussed in a later section. 
The rate, gq, , of production per unit volume of the kth substance may 
depend upon any, all or none of the C;‘ and is negative when the 
substance is consumed. Since q, = 0, for all k, in the environment 
the equation of continuity outside of the cell is 


—V-4°=0, kK=1,2,----8. (12) 


Substitution of J,’ from (10) into (11) gives 
= Di,;' VAG es V Dx;*- V C;' SF a—=—0, k=1,2,--:-8. (13) 
é] 


The second physical approximation which we will make, and which 
is a composite one, is to neglect the second summation in (13). This 
amounts to assuming that the D,;' are: 1) inherently independent of 
the space coordinates and 2) independent of the C;‘. If the Dy;,;' de- 
pend upon the space coordinates only by virtue of their dependency 
upon the C;', the second summation in (13) is 


8 Drji 
a Ci 


For a closed system in which q, = 0, for all k , it is possible to neg- 


p> V Gi V Gf. (14) 
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lect (14) if the initial conditions are that every V/C;' is small, for 
the gradients continually decrease, and if (14) were initially negli- 
gible it may be neglected throughout the course of the diffusion pro- 
cess. This is the familiar proposition of differential diffusion (Har- 
ned, 1947). Of course in such a case the only steady state solutions 
are Cy’ = Cyo' = const. for all k. When every q; is different from 
zero, then, quite regardless of initial conditions, there are steady 
state solutions in which the \/C;' may be by no means small. In fact, 
as will appear in the next section, the vanishing of a particular qx 
does not even insure that the corresponding gradient, VC; , is small. 
Furthermore in most cellular problems the angle between the vectors 
VC,‘ and VC; will be, for all 2 and k, zero or 2. Therefore in neg- 
lecting (14) we must do so on the basis that dD,;;'/oC,' is negligibly 
small and not on the assumption that the scalar products VC,‘ - VC;' 
are small. By assuming, in writing (14), that the D;,;' are functions 
of the C;‘ alone we have also assumed that the interior of the cell is 
a homogeneous structure. This is a virtually necessary, but not par- 
ticularly good, assumption. 

With the above approximations, which are not so readily justi- 
fied here as in some experimental in vitro setups, we take the equa- 
tion of continuity to be 


> Dy; V?2 C5 + q=0, k=1,2,---:8 (15) 
j 
inside of the cell and 
D> Daj? V? CX =0, K=1,2,----8 (16) 
j 


in the environment. Unlike the situation when the simple Fick’s Law 
is employed, the members of the set (15) are not independent, what- 
ever the dependency of the q; upon the C;’. 


Spherical cell; q, constant. 


Consider a cell which is producing s substances at the constant 
rates gq, k =1,2----s. It is assumed that the q are independent. 
Under these conditions of constant q, the set (15) may be reduced 
to s independent equations; each equation is still of second order. 
The solutions of (15) are the solutions of 


V2 C,' =— Q:', 
k=1,2,.-:---8 (17) 
Ox! = Eq; Rei; 
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where R;, is the reduced cofactor of the element Djx' of the determi- 
nant |D,;'|. Similarly the set (16) becomes 


V2 C= 0 pak yas; (18) 
and there are no other solutions, for in general the determinant 


|D;;°| is different from zero. 
The boundary conditions are 


> 


n-dg=de-n,;one; k=1,2----8, (19) 


= 
where o is the cell surface and » is a unit vector with the direction 
of the outward normal to o; 


Peau ie), Paes (20) 


where I, is the flow of the kth substance across the cell membrane 
and 
C= Cen Paste ice, (21) 


For the special case of the spherical cell 
Le (rC;,*) 
CI a 


VG , (22) 


and the solutions of equations (17) are 


f= et Ay, PST, j= 162 ess (23) 


where the A; are constants of integration and 7 is the radius of the 
cell. We have used the condition that C;‘ be finite at r = 0 to elimi- 
nate s of the constants which appear in the integrals of (17). The 
solutions of (18) are 


Cio = Co = Ball, 92 Tepe — AP eee ce (24) 


where the B; are constants of integration and we have used condi- 
tion (21). 


The boundary conditions (19), with (28), (24) and (10), give 


(25) 
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which set determines the B;, as 


Wena To 
B.=—> q; Li; =— Qi? (26) 
Bue 3 


where L;; is the reduced cofactor of the element D jx° In the determi- 
nant |D;;°|; Q:° is the same function of the D,;;° that Q;‘ is of the Dyj'. 

The proof of the last equality in (25) is.as follows: From the 
definition of Q; it follows that 


> Dx;' Qj) = > Di Dj Fein In = > Go One 5 (27) 
j 7 Nn n 
where 6, is the Kronecker delta. The term with n = k is 


OX Di; R= Oe S Dy’ Ryy=Qk, K=1,2----8, (28) 
j j 


because the summations in (28) are the expansions of the determi- 
nant |D;;'| by elements of the kth row or column, divided by the de 
terminant. The terms with n + k are 


dois ig 0, k=1,2,----,#n 
j 


N=1,2,:-::-,#k Bs, 


since the terms in these summations are the products of the reduced 
cofactors of the nth row and the elements of the kth row. By known 
properties of determinants (Aitken, 1948), expansions in terms of 
alien cofactors vanish identically. 

The imposition of conditions (20) requires some discussion: 
The natural and straightforward generalization, which is consistent 
with (10), of the conventional permeability expression requires that 
equations (20) read ) 


0C;' 
or 


—E Dis( ) = Shas (Ci§—Cy) _, =, 2---+8. (80) 
j j 


T=T0 


The right hand-side of equation (30) is I, , the flow across the mem- 
brane in the positive r-direction of the kth substance; the j,k “~7, 
are the “cross permeability coefficients.” The above form of I; does 
not imply any particular mechanism of permeation. It expresses 


the assumption that if the flows, J;' (o), to the membrane are mu- 
tually dependent, then the flows across the membrane also are, and 
in the same manner. In the simple case the flows through the mem- 
brane are proportional to the concentration differences across the 
membrane; in the generalized case the flows across the membrane 
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are linear combinations of the concentration differences. This ap- 
pears to be a necessary assumption, whatever the physical basis of 
permeation, whenever equation (10) is assumed to hold rather than 
equation (4). If it is assumed that the membrane is a distinct “bulk 
phase” then 


(€3-—G;) 


, (31) 


=—Z Dan ) => K; Dx; 
j OFZ] aoa eg 

where the Dyjm are the coefficients in the membrane, (@C;/dr)m the 
gradients in the membrane, K; the partition coefficients between the 
environment and membrane substance and 6 the thickness of the 
membrane. The generalized permeability coefficients are then 


hij = K; Dyj"/6 . (32) 


In (80) and (31) it has been assumed that the partition coefficients, 
a;, between the environment and cell substance are unity. This re- 
striction is easily removed. Needless to state, there is some doubt to 
be raised with regard to considering the membrane [probably 6 < 
10mu (Danielli, 1942)] as a “bulk phase.” It is interesting to note 
that in the “stationary film” theory of the kinetics of gaseous absorp- 
tion into liquid phases, films of the order of 10 cm. (Roughton, 
1941; Davis and Crandall, 1930) have been regarded as “bulk mat- 
ter’ with good agreement between theory and experiment. 
With equations (23), (24) and (26), equations (30) give 


r ‘ ‘. 
D My Aj =— > Diy’ Qi) + ES Ixy D;, 
j 3 Oj j 
wots . io Sie 
Ps eG soles 


Using the result (27), the solutions of equations (33) may be writ- 
ten as 


r 
Aver aie Q; Aig + SD hin Dn Ai; , (34) 
j jon 


where H;; is the reduced cofactor of the element h;, in the determi- 
nant |’.;|. By the same argument, equations (28) and (29), by 
which (27) was established, it follows that 


2 Pn ~ Nin Aix; = 2 Dn One (35) 
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Therefore, with p;, from the definition in (88), equations (34) give 

Yr {Birk To" 
‘A, == (,° + —Qyn ae ey a ae eka k=1,2,-:--d-—1(36) 

3 6 S 

where 
Qn” => a; Ais, b= Then Sees (37) 
wv] 


may be formed from Q;,' by replacing D,;* by hy; . 
If in equations (30) the terms in I; are taken as 


hj (Cj'— 0 C;°) 
2 
then in (36) C;° is replaced by a; C;,° and Hm Q:° is replaced by 
oO 


fee 
Oi eA Ee 
3 


This completes the solution of the problem with the final re- 
sults: 


. r a 2 
Ci Cs + — Q,' 4 a (1.2 — 72) On. 
3 6 3 
| Se ee oc 


(38) 


(Pe 1 
C=C + OQ, b=1,2,---8. (39) 
Oo 


The solutions of the same problem under the assumption of the 
simple Fick’s Law are 


Yo OK vr Yo Qk 

Ci = Co +—— + re Citptensk ARS Meenas 40 

AGA re eye err te See 
Ton Ok 1 

e=C9 4 =, 41 

Ce Oe (41) 


We will refer to this as the simple case as opposed to the generalized 
case (38) and (39). The simple and generalized cases may be com- 


pared as follows: 


1) In the simple case the parameters in C;’ and C;° are deter- 
mined by the properties of the kth solute alone. In the generalized 
case these parameters are determined by the properties of all of the 


solutes present. 
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2) In both cases the C;,' may be separated into terms contrib- 
uted by the membrane, internal diffusion properties and external dif- 
fusion properties. The simple constants hi, Dia’ and Dy? are re- 
placed in the generalized case by the corresponding determinants 
\An;|, |Dxj'| and |Dx;°|; the gq, are replaced by the appropriate deter- 
minants from (17), (26) and (87). The characterization of the 
membrane properties by the set hx; is formally similar to the situa- 
tion discussed by M. F. Morales (1944). In the Morales’ case the 
mutual dependence of the diffusion currents is based upon a physi- 
cal situation somewhat different from that implied by equation (10). 
It is worth mentioning that if equation (82) is accepted, then Q,”, 
from (37), is proportional to the membrane thickness 6. Thus in 
problems (Rashevsky, 1948, p. 108) in which it is necessary to in- 
voke explicitly the dependence of permeability upon 6, the same for- 
mal relations may be preserved in the generalized case. 


3) In the simple case C;'(7”) is a parabola with apex at the 
cell center open downward if q > 0; open upward if gq, < 0: itisa 
constant, (C;° if o,—=1), if q—=0. The differences C;' (0) — Cy! (10) 
= OxV 07 / 6D ix! and C;? (%) — Cx? (7) = 1 o/ dhix » & =—1 , are direct- 
ly proportional to and have the sign of q,. Thus in the simple case 
the concentration of a produced solute always increases toward the 
cell center and is always higher everywhere in the cell than in the 
environment and conversely for a consumed solute; the concentra- 
tion of an inert, q, — 0, solute has a constant value in the cell and 
environment. 

In the generalized case C;'(7”) is a maximum or a minimum at 
ih ee, 0 when Q;' > 0 or when Q,;' < 0 respectively. The differences 
Cx! (0) — Cu! (170) = Qe! 170?/6 and Cyi (70) — Cr® (110) = Qe” 70/8 are lin- 
ear in q, and have the sign of Q;' and Q,” respectively. Thus in this 
case the concentration of a produced solute or of a consumed solute 
may increase or decrease toward the cell center and it may be that 
Ci (70) < Cx? (10) or that Cy' (70) > Cx® (70) ; further, depending upon 
the sign of Q;°, C;,° may increase or decrease toward r = o: Similar 
statements apply to an inert solute (q,= 0). The vanishing of gq; 
does not then insure a zero gradient nor are the signs of the slopes 
of Ci (7) and C,°(r) determined directly by the sign of gq, as in the 
simple case. These points have been discussed elsewhere (Hearon, 
1950b) for the special case s = 2,q,=0,q@ <0. 


When Q;' = 0, C;’ has a constant value in general different 


from Ci even if a, =1. Sucha situation would simulate, under the 
assumption of the simple Fick’s Law, a substance for which a = 0 
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ly Dix! —= ©. Similarly Cie = C,° tf Ox? = 0 and C! (7%) = Ce (ro) 
if Q.” = 0. The difference between the mean cellular concentration, 


1 


Cy ae {| C.'dV, and C;,°, the unperturbed external concentra- 
= 


tion, 
_—. r Y, 2 } 2 
Cy = Co + = Qu” + — Qi + = Q:° (42) 
3 15 3 


is the quantity most likely determined experimentally. This quan- 
tity is positive or negative depending upon the relative values of Q;', 
Q;” and Q,° subject to the restriction on C;° that 

A,>0O, (438) 
which insures C;,' > 0 for 0 < r < 7, which we have assumed 
throughout. It is probable that Q;', Qx° and Q,” will have like signs 
(for the values of the D;;, 7 # k are not unrelated to those of the 
Dx), but this is not a necessity except in special circumstances [e.g. 
the special case already treated (Hearon, 1950b) ]. 


4) The direction of the currents J;' and J;° are in both cases 
determined by the sign of q;, for in both cases the currents are 


Tisstr; (44) 
3 
aoe r 
Jeo = ee (45) 


These relations are obvious for the simple case. For the generalized 
case equation (44) follows from equations (10), (23) and (27); 
equation (45) from equations (10), (24) and (25). In both cases 
then the currents vanish with q,. In the generalized case the gra- 
dients VC,’ and VC;° have directions determined by the signs of 
Q,' and Q,° and do not vanish with q. In fact if the signs of q: 


and Q,' are different, then J;' and VC;' are oppositely directed and 
the solute flows opposite to its own gradient. There are however 
thermodynamic restrictions upon this situation: The quadratic form 


r \2 oe _( os Rij 
( ~) Sac5= c) BE macs 
wd 


= Mi Jess, 
kf 


(46) 


Mij = Ry; /Cx' 
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must be positive definite. For if equations (10) be solved for the 
VC; there results 


V Chi =—D Fi; J;'- (47) 
j 
Observing (44) it is seen that (46) is 
1 dF ) 1 ase 
areal k= ea =————5S V weds 20 (48) 
RT ( dtc)aee RT = YM 


[ef. Hearon, 1950b, equation (35) ], where F is the Gibbs free energy 
per unit volume. The matrix [M;;] is, excepting the scalar multi- 
plier RT , the reciprocal of the matrix [Qj] as is seen by comparing 
(47) with the solution of (1) for Vu and taking u = 4° + RT In C,. 
Similarly, by sclving equations (1) for VC; it is seen that the ma- 
trix [R;;] is the reciprocal of the matrix [0,;/C;‘]RT . Except for 
a scalar factor then, (46) is the dissipation function of Onsager 
(1931a, 1932, 1945). The requirement that this be positive definite 
imposes inequality restrictions upon the D;;' (Onsager, 193la; Ec- 
kart, 1940; Meixner, 1942). 

Equations (17) and (18) and the boundary conditions (19)-(21) 
are independent of the coordinate system appropriate to the prob- 
lem at hand. Thus in the integrals of (17) in any coordinate sys- 
tem (e.g. Young, 1939), obtained under the assumption of (4) and 
gd. constant, one may replace q;/Dyx' everywhere by Q;' and obtain 
the generalized case. The subsequent imposition of (19)-(21) is 
conceivably a difficult problem in general. 

It will sometimes, but not often, happen that certain members 
of the set (17) are integrable by quadratures when one or more 


of tne q, depend upon C;’. We will in fact later discuss one such 
case. 


Spherical Cell; qy = >) kx; C;*. 
j 


When the q are linear combinations of the C;‘, the set (15) 
can be solved, although if s > 2 difficulty of an algebraic nature 
arises, viz., the solution of an algebraic equation of degree > 2s. 
When the q are proportional to a concentration the situation is, of 
course, simpler but even with s = 2 there is considerable algebraic 
tedium associated with the boundary conditions and the results can- 
not be simply and concisely contrasted to the simple case. 

We wish to point out for the simple case s = 2, q, = kG 
Q2 = kzC,' the dependence of the form of the solutions upon the 
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physical situation. The set (15) is written, using (22), as 


> Dyj' D? X; + ky, X,=0, k=1,2, (49) 
j 
where 
Xi, =r Gas 
ped, 2 
dr 


The roots of the auxiliary equation are 


poe a se Woes dil, As 

A 2A aa le) B ’ (51) 
where a = kz Dy' + ky Dz’, A = |Dy;'| and a is associated with the 
positive sign in (51). Now it will be recalled (Rashevsky, 1948) 
that in the corresponding simple case (D,.' = D.,‘ = 0) for a pro- 
duced solute X; = A; sin ar and for a consumed solute X; = A; sinh 
ar. In the more general case X; has a variety of forms for the roots 
(51) may be all real, all pure imaginary, two real and two pure 
imaginary or all complex. For brevity these cases are tabulated be- 
low in Table I giving the form of X; after the elimination, by the 
condition C;' = finite at r = 0, of extraneous constants. It is as- 
sumed throughout that 4 > 0. The constants gi, Qie, Gus are not 
independent. When both substances are produced, or both are con- 
sumed, the form of X; is similar to that in the simple case, pro- 
vided that a? > 4k,k.A, but two discrete sin or sinh terms appear. 
When a? < 4k,k.A , which can occur only if k, and k, have like signs, 
the form 5 holds whether the substances are produced or consumed. 
For the special case a? = 4k,k.A the sin appears if the substances are 
produced; the sinh if they are consumed, but these terms are multi- 
plied by a term linear in 7. If one substance is produced and one 
is consumed, the form 6 holds for both substances. The special case 
a = 0 is presented because of the relation of the roots in this case 
to those in 3 and 4: The form is not distinct from that when a # 0. 
Cases 8 and 9 give the distribution of an inert solute in the diffu- 
sion field of a consumed and produced substance respectively. The 
immediate effect of the inert solute upon the metabolized one is to 
increase 6 since k.D,,'/A > k.D2.' but the values of g. in these cases 
will be different also from those in the corresponding simple case. 


discriminant > 0 
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TABLE I 
CONSTANTS Roots xX 0," 
ke, kg > 0, = 40, Ie SiN AT + Ix Sin BL 
discriminant > 0 ae 4), /6) 
hha <0, Eo Ix, Sinh ar + Jy. Sinh Br 
ae 


levy he: > Os E iy hika/A (Gir + 9x2 7) Sin aT + Oy, 7 COSAaT 
discriminant — 0 HE iy hike/A 
Ke, Ke, <0 Ey haha /A (9x1 + 9x2 7) Sinhar + 9,3 7 coshar 
discriminant = 0 Wy hike/A 
k,,k, $0 a, +716, 9, (e%" sin b, r + e”7 sin byr) 
discriminant < 0 a, = 7b, 
ke S0,k, 20, +a Ix, SNH aT + Jpo Sin BT 
ax0 == ye) 
1. hk, S0,k 20, + y/-haka/K Ix, Sinh aT + Jpo sin Br 
G0 + iy/—hake/B 
3. Ky—= 0, ta —0 X,— 9, sin Pr 
—ke Dy,é A 
k, > 0 + age X,=91,8in Br + 9.7% 
aici see LY oe 0 X,—=g, sinh Br 
—k, D,,* : 
ky <0 + A =+£ X,=—g9,, sinh Br+ 9327 


Coupling by Diffusion. 


It can be seen from the above results that the rates of chemically 
independent reactions are mutually dependent. The rate of a given 
reaction is dependent upon the occurrence and magnitude of the 
rates of other reactions which are proceeding, but with which it 
has no products or reactants in common. Furthermore, from cases 
8 and 9 above it is seen that an inert substance may influence the 
rate of metabolism of a given solute. It has already been seen that 
the effect of the drag coefficients in such a case is to increase fp, 
(cases 8 and 9), but also the constant g2 will be determined in part 
by the unperturbed external concentration C.°. Therefore the rate 
of metabolism of “substance 1” will be influenced by the cross co- 


efficients D,. and D., and the external concentration of the inert 
“substance 2.” 
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A detailed assessment of such effects requires evaluation of the 
constants g;; in cases 1-9 and some concept of the magnitudes of the 
Di; , k #7. We wish to exhibit a simple case in which the phenom- 
enon of coupling by diffusion is particularly transparent from even 
the partially solved boundary value problem. 

Consider the case: s = 2, q, = const., ge =— k.C., ke > 0. In 
these circumstances V/?C,' from (17) is directly integrable and sub- 
stitution of C.‘(r) into V°C,' from (17) gives an equation inte- 
grable by successive quadratures. The results are, after requiring 
that X, — 0, atr=0: 

bes “sinha Git es Dee 0: 


Ci =— -A EE ie : 
; Ey: 6D. Hepes oe oe 


sinhar | Ds‘ qs 


Day Ie’ Cy 


Ky De : . 
eee If go = k. C2*, k2 > O then in (52) and (53) 


where a = 


sinh is replaced by sin. The terms in C,' contributed by substance 
2 vanish with D,.', the “drag of 2 on 1.” When D,,' = D.,' = 0, 
(52) reverts to the simple form (40). The additive contribution in 
(53) from substance 1 is proportional to D.,' and to qi. The addi- 
tive contribution of substance 1 to the total rate of consumption of 
substance 2, 4 a ™? Dz:' q:/3 Di‘, may be appreciable if q: is suffi- 
ciently large. If it is assumed that h.; = 0, which is no essential re- 
striction, the condition (20) becomes for k =2, 


C=C al ei — Fo (54) 


If conditions (19) are replaced by 
[ aav= [ tv-do;  k=1,2, (55) 
Lice oO 

the quantity A, may be evaluated rather directly giving 


Vidi Dy cor 1 | Det a 
NG tr 
a} Dui eS iiawh. - Dulla| 
A. = —— orm? (86) 


° . De Dix? 
sinh a7, + ; a7. cosh a7, — sinh aro eri hs 


where A, = |D;;°|. The constant B, in C2’ is determined from 


GENERALIZATION OF FICK’S LAW 


fat 
Or 
bo 


Bz = As sinh aN SF | ae Coe Toe (57) 


If substance 2 is consumed, sinh and cosh are replaced by sin and 
cos. If gq: < 0, then A, > O and it increases linearly with |qi|; Az 
also increases with D.,‘ but the dependence is not simple since D2,’ 
appears in a and A. 

Thus the rate of consumption, per unit volume, of substance 2 
is influenced by the presence and the rate of consumption of sub- 
stance 1 due to a) the influence of the additive term, b) the depend- 
enc: of As on g, and D2,* and c) the dependence of f upon D;,'. 

When qi < 0 there are restrictions on C.° such that C,* > 0 for 
0 < r < 7; these limit the permissible dependence of A, upon Dy2 
and De, . 

In an entirely similar manner it can be shown that if substance 
1 is consumed at a rate proportional to its concentration and sub- 
stance 2 is produced at a rate proportional to C,', the rate of the 
reaction is determined among other factors by the presence of the 
product (substance 2) and by the external concentration, C.° of the 
product. Thus the concentration of product in the medium influ- 
ences the rate of a kinetically irreversible reaction: i.e. a reaction in- 
finitely far from equilibrium. Consider the set 


2 

> Dj VW? C;* =k, Ci’ ’ 

j=1 

oF (58) 

»> D.;' Si C;* — ky Cc" , 

j=l 
which describes the production of substance 2 from substance 1 by 
a first order irreversible reaction. The solutions are found, from 
(17), to be 


: sinh ar 
Ci =A,——, (59) 
D;,' a Dz,* ' sinh ar 
Ct = A, — aoa) 
2 2 (aie 1 9 (60) 


where a = Vk, (Do.' + Dy‘)/A. From the boundary conditions 


C.' — C.° at r= To» 


—k, Cc dV=n-J,°, (61) 
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the quantity A, is determined as 


0 C,° 
AS ) ee (62) 


: A Doo® + D2? ( 
sinh ar, + ————-|{ a” cosh ar, — sinh ar 
BOTA, Datei \ ae : ses 


the quantity B, from 


A, sinh ar, — CY To— Bit (63) 


As the coefficients D,,', D,.° increase, A, decreases (largely due to 
the increase in a) and the greater the “drag” of the product upon 
the reactant, the slower the rate of reaction: If Dy? = D2,' = Dy»? = 
D.,° = 0, no constants referable to substance 2 appear in A.. 

In summary then, an inert solute, or simultaneous but inde- 
pendent chemical reactions may influence the rate of a given reaction 
and the diffusion properties of the product modify the rate of an 
irreversible reaction. 


Distribution of a catalyst. 


Granting that a given solute is distributed in the cell according 
to the flows and productions of all other solutes, one may ask in par- 
ticular for the distribution of catalytic particles. Such particles are 
conserved in the cellular reactions, and under steady state conditions 
their rate of production is zero (provided no synthesis of new par- 
ticles is occurring when no steady state is possible). Their dis- 


tribution is determined by the currents, J;', of all other solutes; 
these currents are in turn conditioned by the metabolic rates which 
in turn are determined by the concentration of catalyst. With sim- 
plifying assumptions this problem fails in biological importance for 
the intuitive supposition is verified that if the “drag of outflowing 
solutes” predominates, the catalyst is concentrated at the periphery 
of the cell and conversely. 
Consider for example the set 


2 
“ Dz! Vv? C;' = 
j=l 


(64) 
>» D.;' Wie C;* a he C;', 
j=l 


which describes a substance being produced at a rate proportional 
to the concentration of catalyst (substance 1). From (17), V? C;' is 
integrable giving 
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sinh ar 
r 
Cs Oars ae (65) 
sin ar 
tf 


where a = VDy"|k:|/4 and the sinh form holds for k, > 0. Substi- 
tution of C,‘ from (65) into V?C,' from (17) gives 


sinh ar 
Dui r 


C,§ = A, — Ai ; , 
Dy» sin avr 


(66) 


T. 


The distribution C,‘(7) of the catalyst simulates that of a substance 
which obeys the simple Fick’s Law, is produced at a rate k. Ci‘ and 
exhibits a diffusion coefficient 4/D,.‘. Either of the forms (65) is 
also observed if several substances are involved. For example if 
q:=0, g2=— ke Ci’, G2 = ke C,', ky > 0, i.e. substance 2 is converted 
to substance 3 at a rate proportional to the concentration of cata- 
lyst, then 

WV C= CF el ha els (67) 
and r C,‘ is of the form sin ar or sinh ar depending upon the signs 
and magnitudes of the reduced cofactors R,., R:;. The catalyst is 
not necessarily concentrated toward the periphery or center when 
there is a net outflow or inflow of material 


(2 J, aaree) 


respectively: The distribution is determined by the relative values 
of the D;;' of the outflowing and inflowing solutes. In the case de- 
scribed by (67) there is clearly no net outflow of material for 


5 { a. dV =0. 

k V 
It is not difficult to show that for s solutes with gq, =—0,k=1 nee 
-+++ 8/3 for every catalyst and qx = — dra = kx Ci’, k > 8/3, for 


every product-precursor pair, the conclusion from (67) is general. 
It will be obvious that whatever the spatial distribution of the cata- 
lyst, the total rate of the catalyzed reaction is the same. It is in fact 


kn where n is the total number of moles of catalyst, i) Cav = 1 
4 
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The A, of (65) or the constant of integration in the solution of (67) 
is determined by this condition of conservation of catalyst. This con- 
clusion regarding the equality of total rate for any catalyst distribu- 
tion and conclusions from (64) and (67) [that Cy'(r) = (sinh ar) /r 
or (sin ar)/r] are not true if the rates are not proportional to the 
catalyst concentration, e.g. if gq = k(C,')2, q; = k Cy’ C;*, ete. A 
physically simple example is 


S D,;* ye Cc; =0, 


j=l 


(68) 
= LI 2 CP — ies Cs Gat 
The first of these sdiditte requires that 
D,,' Cy’ + Dp‘ C= A = const. (69) 
and using (69) and (17) gives, 
V2G,' = KC — KC,’ Ce, (70) 


where K, = k, A/A, Ko = kz D,,'/A. If n is smallkthen C,' obeys 
approximately 
d?(r C,') 


= 2 4 Ks 2 
dr —Ki7 Cy $ (7) ? : (71) 


where $(7r) = sinh ar or sin ar if k, > Oor ko < 0,a= \V/|Ki| and 
K,; is K, times an arbitrary constant. The solution of (71) cannot 
be obtained in closed form. An approximate solution is, for k, > 0, 


r Cy =A, sinh ar f(r parece A817) : 
na (—a)”™ orn 
1 ai + (8a)/—2(—a)! 
Cacia 
while the difficulties of the diffusion equation forbid a concise assess- 


ment it is clear that in physically realistic cases a catalyst may ex- 
hibit a distribution different from A (sinh ar)/r or A (sin ar) /r and 


2 (72) 
f(r)= Zar, a= 


that the total rate, ir q dV, of the catalyzed reaction wall depend up- 


on the spatial distribution of the catalyst. The sere of cou- 
pling by diffusion is greatly enhanced by the possibility of distribut- 
ing catalysts in the resultant diffusion field of reactants,. eR EO Ga 'S 
and catalysts. 
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Non-Isothermal Systems. 


The equations of transport for non-isothermal systems are 
available (Eckart, 1940; Leaf, 1946). A few remarks relevant to 
metabolizing systems will be made. It is probable that the gradients 
in temperature due to gradients in composition (inverse Soret ef- 
fect) are small and that the effects of thermal diffusion are likewise 
negligible under these conditions. However a system in which chemi- 
cal reactions are occurring must be treated as having continuously 
distributed sources and sinks of heat. The temperature gradients 
due to this effect are probably important. This effect may be included 
in the problem as follows: In equations (10) » (15) and (16) let the 


summations run from 7 = 0 to 7 = s and let J, be the vector of heat 
flow. Then Dy will be the coefficients of thermal diffusion, Do. the 
coefficient of thermal conductivity and D., the “inverse Soret coeffi- 
cients.” For qo, the production per unit volume of heat, we have 


G=DAa, 
k= 


where H; is the partial molal enthalpy of the jth solute. By express- 
ing the q; from the usual kinetic expressions gq) may be written 
(Hearon, 1950a, b) as 
qo = >> AH i Uj, 
j=1 

where AH; is the heat of reaction and v; the velocity for the jtn 
chemical reaction. Similarly, in equations (30), hro and ho, are the 
coefficients of thermal diffusion and the inverse Soret coefficients in 
the membrane and hy». the heat conductivity of the membrane. Strict- 
ly speaking the q; and v; are, through the temperature, functions of 
the space coordinates and similarly for the 4H; which are tempera- 
ture dependent. This complicates matters considerably but some of 
the problems, involved in simple cases have been solved by W. S. Hor- 
ton (1948). It might be predicted that if q > 0, including heat 
flow in the problem will have about the same effect as an additional 
solute produced at the rate gq. and exhibiting cross coefficients Dox , 
Dio, k =1,2,-+--. 8, There would be superimposed upon the effect 
of coupling by diffusion the effect of “coupling by thermal diffu- 
sion.” It may then be possible to attribute special roles, of mainte- 
nance of non-equilibrium distributions and of coupling, to reactions 
with conspicuously large AH’s. A further implication is that the 
temperature dependence of reaction rates in the intact cell would be 
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a composite affair including the effects of coupling by diffusion and 
thermal diffusion as well as the innate temperature dependence of 
rate processes. Also to the extent that such effects are measurably 
important the frequently determined increase in heat production 
above the “resting level” or “non-nutrient level” will be a complex 
quantity. 


Discussion. 


It may be briefly pointed out that few cases can be solved when 
the gq, are functions of several C;’. When the q; are linear combina- 
tions, p> ky; C;*, many such problems may be solved. Such a linear 


I 
combination suffices, for example, to describe a series of first order 
consecutive reactions and it is easily shown that the determinantal 
equation to be solved is always reducible to degree 2 (s—1). 

The expressions used for the q; in this paper for discussing the 
distribution of a catalyst are most inadequate and the situations dis- 
cussed highly oversimplified. In general the rate of a catalyzed re- 
action is proportional neither to the concentration of catalyst nor to 
the product of the concentrations of catalyst and reactant (Hearon, 
1949). In semple cases the above are limiting cases for very high or 
very low concentration of reactant: In general this is only approxi- 
mately true (Hearon, 1949). 

An important consequence of Onsager’s Generalization is that 
inert solutes may have non-zero gradients and different cellular and 
environmental concentrations. Thus inert solutes can contribute to 
osmotic forces and to diffusion drag forces [if the magnitude of 
these forces be taken as proportional to the gradient (Rashevsky, 
1948; Landahl, 1942)]. Quite different expressions for the force 
components in “Betti’s Theorem” treatment of cell elongation 
(Young, 1939; Rashevsky, 1948) result from these considerations. 

It has been found possible for metabolized solutes to flow against 
their own gradients and for inert solutes to flow into a region of 
higher concentration from one of lower concentration (cf. discussion 
Hearon, 1950b). We wish to point out here the similarity between 
these situations, direct consequences of equation (10), and those 
cases treated by J. Frank and J. Mayer (1947, especially their dis- 
cussion pp. 301-302) and the “diffusion retrograde” discussed by 
T. Rosenberg (1948). 

Finally, in addition to the implications already mentioned which 
equation (10) has for the permeability concept, it would seem that 
attempts to compute absolute permeabilities (Bloch, 1944) must pro- 
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ceed along fairly different lines to the extent that the cross coeffi- 
cients are important. 
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Theoretical considerations lead to the calculation of the amounts 
of air-borne materials removed by the various regions of the human 
nasal passages. These calculations are presented and discussed in the 
light of available data and lead to the conclusion that impaction due to 
inertia is the principle factor in the removal of the air-borne particles 
of diameters greater than one micron. 


In considering the mechanisms by which particles may be re- 
moved from the air stream as they pass through the nasal passages 
(cf. Proetz, 1941; Davies, 1946), one is led to the evaluation of the 
following: impaction on the inner walls, impaction on the nasal 
hairs, settling and diffusion. For all flow rates, and for particles 
greater than one-tenth of a micron, calculations show that the effect 
of diffusion is negligible compared with the effect of settling. 

Schematic Representation of the Nasal Passages. In order to 
make the calculations we adopt the following schematic description 
for the nasal passages. Each of the external nares (A) is taken to 
have a cross-sectional area of 0.75 cm.?, with hairs of about 100 
micron diameter occupying one-half of the projected area. The pro- 
jected area of an object is the area of its projection on a plane per- 
pendicular to the direction of air flow. The projected area of the 
hairs was estimated in one subject by comparing the blackness of 
discs of photographic paper placed behind the hairs. Another esti- 
mate was made by removing all the nasal hairs possible and meas- 
uring the lengths (total length ca. 50 cm.) and diameters. The hairs 
taper and vary greatly in diameter. The variability can be taken 
into account to a first approximation by considering the hairs to be 
of two sizes, one of 70 and one of 130 microns in diameter, each 
group occupying one-fourth of the projected area. 


*The work described in this paper was done in part under contract between 
the Medical Division, Chemical Corps, U. S. Army, and The University of Chi- 
cago Toxicity Laboratory. Under terms of the contract the Chemical Corps 
neither restricts nor is responsible for the opinions or conclusions of the author. 
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The second region (B) is the constriction about 2 cm. behind 
the opening of the external nares at the end of the nasal bones. This 
we take to be a rectangular tube 1.2 cm. high and .25 cm. wide which 
bends at an angle (6;) of 30°. The length of the region is about 1 
cm. 

The third region (C) is taken to be equivalent to a rectangular 
tube 2mm. wide, 3 cm. high and about 1 cm. long, in which there 
is a bend of about 20° (@¢). 

The fourth region (D) will be sub-divided into two regions 
which are connected in parallel, one representing the somewhat nar- 
rower (1 mm.) and more tortuous upper passages, the other the 
wider (2 mm.) and more direct lower passages. Thus region D is 
represented by a rectangular tube whose width is from 1 to 2 mm., 
height about 4 cm., and length about 5 cm., so that the total effec- 
tive wall area for both nasal passages is about 80 cm.” Thus we as- 
sume, on the average, nearly two-thirds of the total area (Hellmann, 
1927) is functional at any particular time. In the upper region some 
parts of the air stream may make several fairly sharp bends. In the 
lower region some parts of the air stream make only a slight bend. 
We shall consider the upper air stream to make the equivalent of one 
change in direction, the angle being about 45°, while the lower part 
makes an average turn of about 20°. Furthermore, we shall con- 
sider our schematized region D to be orientated at an angle of about 
45° from the vertical (y = 45°). We shall make the calculations for 
the case in which half of the flow is through passages 1 mm. and 
half through passages of 2 mm., the proportions of the height per- 
taining to each being in the ratio 3 to 1. The results are rather in- 
sensitive to the value of this ratio. 

The order of magnitude of the pressure drop for the model is 


the same as that observed in the subjects used in the experiments to 
be discussed. 


; Removal by Impaction on Nasal Hairs. The fraction of par- 
ticles moving toward a cylindrical object which is removed by im- 
pingement depends upon the inertia and sizes of the particles and the 
size of the object. We introduce the following notations: The quan- 
tities f’ and f” are the fractions of the projected area covered by the 
two groups of nasal hairs (one-fourth for each group), p and d are 
the particle density and diameter (cm.), V is the average velocity 
of the air stream (cm./sec.) and D’ and D” are the diameters (cm.) 
of the two groups of hairs. Then the fraction of particles entering 
the nose which impacts against the nasal hairs can be written as the 
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sum of two similar terms P;’ and P,,”, one for each group of hairs. 
Thus we write the equation (Landahl and Herrmann, 1940) 
Py =f 150 (p d?V/D')# es | 
O70 aa (pe VLD)? 150'(507V/D7)*_.D’ 

The expression for P;,” is obtained from expression (1) by re- 
placing f’ and D’ by f” and D". To obtain V, the flow rate F' (cc./sec.) 
through both nostrils is divided by twice the area of each nostril. 
Since the velocity distribution is actually not constant, we shall take 
this into account, to some extent, by increasing the variability in 
hair diameter. Thus in the calculations the D’s are taken to be 50 
and 150 microns respectively. The total probabilty, P,, if the hairs 
were regularly spaced, would be the sum of P;’ and P,” for the two 
groups, since these processes are mutually exclusive. To take into 
account the irregular spatial distribution of the nasal hairs, we shall 
introduce P, —1— e’?* as the net probability of removal in region A. 


(1) 


Removal by Impaction Within the Passages. For the probability 
of impaction in bent tubes with parallel walls, and in the absence 
of turbulence, we may derive an expression in the following way. 
Let @ be the angle of the bend. Consider a particle at a distance y 
from the center of the stream approaching the bend. In order to 
reach the outer wall, the particle must slip a distance (W/2 — y) 
cosec 6 perpendicular to the direction of the air stream. A more de- 
tailed calculation for 6 = 2/2 shows that the slip is only slightly 
greater than that which would occur if the particle was shot with 
the same velocity into quiet air (cf. Landahl and Herrmann, 1950, 
Appendix). We shall ignore this discrepancy. By Stoke’s Law, the 
distance that a particle travels through quiet air is p d? Vy/18 » 
(7 =1.8 X 10 poise) when its initial velocity is V,. Since the ve- 
locity distribution is parabolic and the velocity is zero at the walls, 
V, = 1.5V (1 — 4y?/W?). If this distance of slip perpendicular to 
the air stream for a particular position y = y* is set equal to 
(W/2— y*) cosec 6, then the resulting equation determines this par- 
ticular y value so that no particle which has a position y < y* before 
reaching the bend can reach the outer surface. If we solve for y* we 
obtain 

SaaS Welt op 

4 "3 pd? V sin é 
The fraction I’ of particles removed is then obtained by integrating 
V, dy from y* to W/2 and dividing by the same integral but with the 


(2) 
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limits —W/2 and W/2. If pd? V sin @ > 3  W, we find the following 
expression for the fraction of the air-borne particles removed at a 
bend when the flow is smooth: 


2 3 

pie ae eeeae 
12 \ pad’ V sin? 108 \ pd’? V sin 0 

If pd? V sind < 37 W, then —0O. 

On the other hand, if all the particles moved with uniform ve- 
locity, they would be displaced uniformly perpendicular to the air 
stream as they passed the bend, and the per cent arriving at the 
outer surface would be proportional to p d? V sin 6/W. When 
p d@2 V sin 6/187 = W, the innermost particles would have just 
reached the outer wall. Hence the fraction J” removed under these 
conditions is given by 


oie pd Ve sine pd’ V sin 6 


— 4 
18 »W ‘ W (4) 


unless the quantity on the right is greater than unity, in which case 
I" is equal to unity. In regions B, C and D the Reynold’s numbers 
are about 10° to 10? for flow rates under consideration. Since the 
lengths are relatively short and the passage walls are uneven, the 
flow cannot be expected to be smooth. We shall, therefore, use the 
average of the values obtained from equations (8) and (4). 


Removal by Sedimentation. This effect will be considered only 
for the region D. It is appreciable only at small flow rates. For the 
calculations we shall use equation (4) of a previous paper (Landahl, 
1950), with the following change. If every air-borne particle in a 
horizontal cylinder fell vertically a distance equal to the cylinder 
radius, 61% of the particles would reach the bottom surface; in a 
slot the corresponding percentage would be 50. Thus instead of 2.3 
cos y we shall use 2.3 X (0.50/061) cos y = 1.9 cos y. Let 7 be the 
mean time for the air to pass through region D. Then, except for 
particles of diameter considerably less than 1 micron, we may write 
as follows for the fraction settling in region D , 


S i 1 — e73.8x105pa2rcosyy/W (5) 


Let A; be the surface area of region D. The quantity F is the 
flow through both sides of the turbinates, W is an average width of 
this region. Since F = 2H W V,7 = L/V and A; = 4HL, where H 
and L are the height and length of the region, the time of passage 
can be written as 
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7 AGW OF. (6) 


Introducing this value for +r in equation (5), we find, since angle 
y= 45°, 
Sa = 1 — E77 -6x108A p@y/F (7) 


Calculation of Total Removal. For region A the net probability 
of an air-borne particle being removed is Py = 1 — e-’*, where P, = 
P;, + P,, the sum of the probabilities for the two groups of nasal 
hairs, if we neglect impaction in this region. For regions B and C, 
P; =I, and Py = I,, since only impaction is considered in these re- 
gions. For region D we have Pp = Ip + Sp — IpSp. Thus the total 
probability of removal can be written: 


The amounts removed in the various regions may also be obtained. 
TABLE I 


CALCULATED PERCENTAGES OF MATERIAL REMOVED IN VARIOUS PARTS 
OF THE NASAL PASSAGES FOR DIFFERENT PARTICLE SIZES 


Diameter 
microns RegionA  RegionB RegionC Region D Total Total 
(Imp. Eq.) (Nasal Hairs) (c=) (e179) 
(4.5 1/min.) 
2 7 0 0 1 5) 6 
5 9 1 1/2 3 14 20 
10 25 3 1 10 39 46 
20 40 22 2 21 86 89 
30 47 40 4 8 99 99 
40 52 47 1 1/2 100 100 
(18 1/min.) 
1 1 0 0 0 1 2 
vs 4 if 0 1 6 8 
3 11 i 1 2 14 17 
5 23 3 il 3 30 34 
7 31 5 2 10 48 53 
10 36 23 3 16 78 81 
15 42 44 5 a 98 99 
20 45 53 1 0 99.9 99.9 
25 48 52 0 0 100 100 
(72 1/min.) 
1 3 1 0 1 5 6 
2 17 te 1 2 ie 23 
5 39 21 2 14 TEL 78 
10 41 57 1 1 99:9 99.9 
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In region A the proportion is just Ps. In B it is (1 — P4)Ps; in 
CG it is (1 — P.) (1 — Ps) Po; and in D itis ((1 P,) (1 — Ps) 
(1 — Pco)P». The sum of these can be checked by equation (8). 

Since the density and particle diameter occur in the combination 
p d? with but one exception, one would expect it to be convenient to 
introduce the quantity d’ = \ pd? in both the calculations and the 
handling of data. This quantity has been referred to as an equiva- 
lent impaction diameter (Landahl and Tracewell, 1949). This quan- 
tity will be used throughout. 

Table I gives calculated values for the amounts of various sized 
particles removed in various regions at three flow rates. These val- 
ues give some idea of the shift in locus of deposition which can be 
expected with change in the flow rate. Also, in the last column are 
values for particles of density of 1/9; these give some idea of the 
expected effect of change in density, e.g. light flocculent material. 

The results of these calculations may be compared with avail- 
able data (cf. Landahl and Tracewell, 1949). We consider first the 
effect of the nasal hairs. The comparison is given in Table II. 


TABLE II 
REMOVAL OF PARTICLES BY NASAL HAIRS 
(Flow Rate—18 l/min.) 


Cloud Fraction Diameter Calculated Removal Observed Removal 
(Imp. Eq.) % Jo 
A 13 43 (78) 
B 8.7 39 44 
C 3.9 27 24 
D 2.0 9 2 
Imp. 1.1 1 (8) 


The comparison for the total removed at 18 l/min. is shown in 
Figure 1. It will be noted that the principal deviation is that the cal- 
culated values of the per cent penetrating the nose are high for par- 
ticles of about three microns in diameter. If we go to higher flow 
rates (cf. Landahl and Black, 1947) for corn oil, we see that the 
shift to smaller sizes is somewhat higher than expected—abosut two 
and one-half instead of almost two. Since for flow rates appreciably 
greater than 18 l/min. sedimentation is negligible, as is the factor 
d/D in equation (1), the effect of a change in flow rate by a factor 


of two should be equivalent to a change by a factor of \/2 in the par- 
ticle size. It should be pointed out, however, that the reflex changes 
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EQUIVALENT !MPAGTION DIAMETER IN MICRONS 


FIGURE 1. Per cent nasal penetration of various air-borne particulates: 
O—corn oil, N—NaHCO, , C—Ca,(PO,),, M—methylene blue, B—bismuth sub- 
carbonate, T—tyrosine, S—Na,SO,, G—glycerol (Landahl and Tracewell, 
1949). Broken line: theoretical curve from calculated values in Table I. 


which occur upon marked changes in flow rate, temperature or hu- 
midity make it difficult to make accurate comparisons. 

The effect of subject variability and variability among subjects 
would tend to lower the high values and raise the low values where 
the curvature is great. An estimate of this effect suggests that it is 
appreciable. Furthermore, the method used (Landahl and Black, 
1947) in which cloud fractions are divided, does itself also introduce 
variability, so that actually for a single subject over a short time, 
the shape of the curve is more like the shape of the calculated curve. 
In the case of non-liquid materials, the failure of the particles to ad- 
here perfectly in the sampling instruments adds further variability. 
Also the failure of some of the particles to adhere to the hairs or the 
walls of the nasal passages, as well as the possible breaking off of 
parts of aggregated material, tends to increase the variability. It 
will be noticed that the corn oil gives the sharpest curve. The points 
representing glycerol should not be considered. In this case the ef- 
fect of humidification on such a hygoscopic material cannot be neg- 
lected since the particle may increase in size due to absorption of 
water. Humidification may also play some role in the case of some 
of the other substances. 
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We mention, in passing, some factors which have not been con- 
sidered but which could conceivably influence the percentage re- 
moved. If the nasal linings are considered to be sufficiently good 
conductors, the only way in which the electrical charges on particles 
could operate to influence removal would be through repulsion be- 
tween like charges. However, it is likely that, unless a special effort 
were to be made by the experimenter to reduce the time between for- 
mation and sampling of the charged particles, the effect of such re- 
pulsion would be spent before the particles reached the nostrils, since 
the time of passage is a small fraction of a second. 

The effect of the temperature gradient in the experiment dis- 
cussed is probably very small since the actual temperature difference 
between the incoming and outgoing air was seldom much more than 
10°C. 

The effect of humidification on the removal of non-hygroscopic 
particles can be roughly estimated in the following way. For sim- 
plicity consider a right angle bend and take a particle at the center 
of the stream which would, due to impaction, just reach the outer 
wall. If the humidification were taking place uniformly in region D , 
then the velocity uw of flow perpendicular to the air stream would be 
u=qF/4H L, q being the fraction of the volume occupied by water 
vapor. For q we shall take 0.03. Although the actual value is gen- 
erally somewhat larger, the viscosity of water vapor is smaller than 
that of air. We shall thus not introduce a change in the viscosity. 
Instead of a movement by the particle into still air we shall now have 
a movement into air which is moving in the opposite direction. Thus 
we have, if the mass of the particle is m, 


mz=—B8and(4+u), (9) 
so that for = 0 and x =1.5V at t = 0, we obtain 


LbV + 4 = 
e=—ut + Sata m( 1 e€ Ne (10) 
Ifu=0 and t= o, then « = x = p(3V/2)d?/18y , the expression 
used in obtaining equations (3) and (4). 

The maximum distance x,, of slip through the air stream occurs 
at t* when the velocity is zero. If the wall is not reached, the par- 
ticle then begins to turn back because of the velocity u. Setting the 
velocity equal to zero, solving for t* and substituting the result into 
(10) we find the maximum distance of slip to be 
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tn = 2 ( aa einen 11 
miata acs 10g | |. 

3V ~ 2u ) ue 
Since F = 2V WH, we find the following expression for the relative 
error made in neglecting this effect: 


(Xo — Lm) /Lo = —— log —. (12) 


If now q = 0.03, L=5, W = 0.15, then the relative error is 0.002. 
The relative change in the percentage removal is thus negligible in 
this case. Where the percentage removal is small, the relative error 
becomes large, but the net effect is still small. 

In the light of the above discussion and in view of the crudeness 
of the model, the calculated results are fairly satisfactory and are a 
further demonstration that the principal factor in the removal of 
air-borne particles by the nasal passages is impaction due to the in- 
ertia of the particles. 

The author wishes to express his appreciation to Dr. H. B. Perl- 
man, Associate Professor of Otolaryngology, for reading and discuss- 
ing the manuscript. 
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If an enzyme acts in a manner somewhat different from that gen- 
erally considered, it may appear to act as a pair of enzymes. This 
similarity may remain after various changes in the solution such as 
activation by magnesium ions. Conditions are given under which this 
type of activity may be differentiated from that of an enzyme pair. 
Some alkaline phosphatases exhibit the type of time course of hydro- 
lysis under consideration. Available data do not eliminate the possi- 
bility that some alkaline phosphatases act in the manner here suggested. 


The study of some alkaline phosphatases suggests that two 
closely associated, inseparable, enzymes are present in these cases 
(Householder and Gomori, 1943). The purpose of this note is to 
consider some conditions under which a single enzyme may give re- 
sults similar to those obtained from such an enzyme pair. We thus 
consider the following scheme (H = enzyme, OP = organic phos- 
phate; cf. Schwab, 1937, p. 337) 


ky Ke 
Per OR OPE. OL +P; 
Kx 
ks 
Of > OnE x. 


We shall consider the last two steps to be irreversible since in this 
system the addition of products does not have an appreciable effect 
up to one hundred times the concentrations normally used. Further- 
more the concentration of OP is assumed to be constant because only 
very little material is broken down. 

From the above scheme we may write, if S= [OP], x = [OPE], 
y = [OE] are the corresponding concentrations, 


oe ES — kaw — kat; (1) 
Yo hege — Tey; (2) 
dt 
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and the total enzyme concentration is 
Eoe=2 yes. (3) 


If we denote by é and 7 the deviations of # and y from their steady 
state values, we have 


&=a4—kh,k,ES/[kikS + kiksS + kaks + keks] (4) 
and 
n= Yy — kik. S/[kikeS + kik3S + kaks + keke] ; (5) 
then 
dé 
pee ee (KS + ka + ke) & — kSn; (6) 
haf — han. (7) 
Thus the solutions for € and y are of the form 
E=Ae™' + Be**; (8) 
n= Ce + De; (9) 


where 


Ax,2 =4(k,S =e hese Se ke 45 ks) 


So Se ee ee 
+ $V (k,S She kK, ta Ky AF a2 ti Cie hee Ee k,Sks a5 kikeS =e Kknke). 


Experiments are generally carried out with rather large 
amounts of the substrate. Thus if S is very large one of the com- 
posite rate constants, 4,, becomes infinite; the other, 4, becomes 
k, + ks. If also, y= 0 and x = E, for t =0, then we find from (3), 
(4), (6) and (8) that 


k3Ho kK o 
— Jeb pe cemneanecg 
kes te ks ke ae ks 
Since the rate of production of inorganic phosphate is kx, we have 


for the total amount of inorganic phosphate [P] at the time t, if 
ak, + k, andif [P] =0fort=0, 


es Kaleskict k?Eo 
So a” 


e7 (katks) t zs (11) 


[P] Oe td fe (12) 


If we set 
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kB 
Ur 2 5 
k. ames i (13) 
ee 14 
Ur 3 ke a ke ’ ( ) 


then equation (12) becomes identical with equation (6) of A. S. 
Householder and G. Gomori. Since generally it was found that 
XL; >> y,, this would mean that k. >> k,;. Thusa =~ k.. 

Where the effect of storage (Householder and Gomori, 1943) is 
such as to affect both components in about the same manner, the in- 
terpretation in terms of the present scheme is evident. In cases in 
which this is not so, it is necessary to suppose that incubation has 
caused a differential change in the reaction rates. 

The effect of magnesium ions on the activity of certain enzyme 
preparations at low values of the pH (ca. 8) can be accounted for 
by considering that addition of the ions greatly increases the rate 
constant k; but only slightly increases k.: thus the “x” component is 
apparently inactivated while the “y” component appears to be acti- 
vated. At higher values of the pH (ca. 10), both rate constants are 
slightly increased by magnesium ions. 

When the enzyme preparation is incubated for several days with 
the substrate containing magnesium ions, the activity generally ap- 
pears to fall off slowly instead of maintaining a constant activity. In 
terms of the two-enzyme hypothesis, this would be interpreted as 
meaning that the “y’’ component is very slowly inactivated. In the 
present scheme it would mean that either of the enzyme complexes 
containing (or in the presence of) magnesium is inactivated very 
slowly. Let the inactivation rate constant for the (OPE) complex 
be m* and let the rate constants in the presence of magnesium ions 
be denoted by asterisks. In this case we find an expression for the 
amount (P*) of inorganic phosphate formed in the presence of mag- 
nesium ions in a manner similar to that used to obtain equation (12). 
Since m* is very small, we may simplify the expression, in which case 


we have: 
ke (1 as 2m*k*; ) ( 1— ene ) “5 ies ( 1 
m* a* at? 


eee | x ( L— en (atem*h*2/a*) t ) : 
a*? 


Pas 


(15) 


This expression is entirely similar to that which is obtained when 
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the “y’’ component is taken to be inactivated [Householder & Gomori, 
1943, equation (5)]. However, the interpretation of the data may 
be simpler in terms of one of the hypotheses. 

We next consider the possibility of differentiating between the 
two hypotheses by experiments involving addition of a tagged sub-. 
strate at various times 7’ with subsequent determination of the rate 
of hydrolysis of this substrate at the time of its addition. Suppose 
that k, and k., are much greater than any other rate constant. Then 
if the new substrate S’ is added in such an amount that after addi- 
tion the concentrations of S’ and S are in the proportion S’ = RS, 
then we will have «’ = Rx, the primes denoting the tagged mole- 
cules. Since {d[P’]/dt},, = k. a’, the rate is also equal to kz Raz 
and thus {d[P’]/dt},, = R {d[P]/dt},-. 

Experiments show that phenylphosphate and f-glycerophos- 
phate appear to be almost interchangeable for one enzyme prepara- 
tion. Thus suppose that the enzyme preparation has been incubated, 
with and without magnesium ions, for a time 7’ with $-glycerophos- 
phate as substrate. Then the addition of phenylphosphate at time 7’ 
allows one to determine the amount of phenol produced within a 
short time after 7’ and thus {d[P’]/dt}, is calculated. But since this 
is equal to R{d[P]/dt},. , this measured rate is just a measure of the 
slope of the original curve at time 7’. If k, and &, are not suffici- 
ently rapid this interpretation cannot be used. The ratios of the 
rates obtained with and without magnesium are the magnesium ac- 
tivation ratios for the corresponding times 7’. Experiments of this 
kind (Gomori, personal communication) show that the rate 
{d[P’]/dt},, is large for small 7’ and it decreases to a small constant 
value in about the expected manner. In the presence of magnesium 
ions the rate is about fifty per cent larger initially and decreases 
thirty per cent in twenty-four hours. The magnesium activation was 
about forty-four per cent for the first hour and slightly over seven 
fold at twenty-four hours. This can be accounted for in the present 
hypothesis by setting k*./k, = 1.1 and k*,/k; = 10, since a = 0.9 
Hr and. he / k= 20 

It can be seen from the reaction scheme that, if the substrate is 
removed, “x” and “y” are eventually converted back to E. To in- 
sure that this will take place, the solution of enzyme and substrate 
should be dialyzed, for a time much greater than 1/k;, against the 
same buffer used for the hydrolysis and without change in tempera- 
ture. This time, 1/k; , is roughly equal to (#-/y;) (1/a). If the time 
of dialysis is +”, then upon the addition of substrate equation (12) 
again holds, but with «, replaced by 2-(1 — e*7”), Since 1/oa is not 
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more than a few hours while x-/y; is not likely to be much more than 
one hundred, the original situation should be restored after a few 
weeks. 

In a preliminary experiment, which was part of a program to 
test the hypothesis presented here (Gomori, personal communication) , 
enzyme preparations were allowed to act upon the substrate for a 
time much greater than at after which they were dialyzed up to four 
weeks against distilled water. The addition of buffer and substrate 
resulted in a curve for the hydrolysis which, after a twenty minute 
delay, showed a slight inflection and then continued to rise at a con- 
stant rate typical of the “y’’? component. However until the results of 
experiments are available, in which the dialysis is against the buf- 
fer, it is not possible to exclude the possibility that, in the case of 
this type of alkaline phosphatase, a single enzyme is acting in a man- 
ner, as here suggested, so as to appear to be two components. 
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